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ABSTRACT 

Wavelets on the sphere are reintroduced and further developed independently of the original group theoretic 
formalism, in an equivalent, but more straightforward approach. These developments are motivated by the 
interest of the scale-space analysis of the cosmic microwave background (CMB) anisotropies on the sky. A 
new, self-consistent, and practical approach to the wavelet filtering on the sphere is developed. It is also 
established that the inverse stereographic projection of a wavelet on the plane (i.e. Euclidean wavelet) leads to 
a wavelet on the sphere (i.e. spherical wavelet). This new correspondence principle simplifies the construction 
of wavelets on the sphere and allows to transfer onto the sphere properties of wavelets on the plane. In that 
regard, we define and develop the notions of directionality and steerability of filters on the sphere. In the 
context of the CMB analysis, these notions are important tools for the identification of local directional features 
in the wavelet coefficients of the signal, and for their interpretation as possible signatures of non-gaussianity, 
statistical anisotropy, or foreground emission. But the generic results exposed may find numerous applications 
beyond cosmology and astrophysics. 

Subject headings: cosmology: cosmic microwave background — methods: data analysis 



1. INTRODUCTION 

The last decade has recognized the cosmic microwave 
background (CMB) as a unique laboratory for achieving 
precision cosmology. The cosmological parameters defin- 
ing the structure, the energy content, and the evolution of 
the un iverse are now determined with an impressive pre - 
cision (Page et al. 20031 ISpergel et al. 2003t iBouchet 2004). 
However, the theoretical hypotheses on which the corre- 
sponding concordance cosmological model relies still must 
be fully investigated. They notably extend from the cos- 
mological principle of homogeneity and isotropy, to the 
model s of inflation for the physi cs of the early universe 
(Wand elt2004 iBartolo et al. 2004 . o r the theory of gravi- 
tation itself, namely g eneral relativity JBoucher et al. 2004a: 
Boucher et al. 2004-3). 

In the context of the concordance cosmological model, the 
cosmic background radiation is understood as a unique re- 
alization of a gaussian and stationary random signal on the 
sphere, arising from quantum energy density perturbations de- 
veloped in a primordial inflationary era of the universe. In 
this respect, the analysis of the CMB anisotropies on the sky 
is essentially confined to the study of its temperature (and po- 
larization) angular power spectrum. 

But questioning the basic hypotheses of inflation, or of the 
cosmological principle notably amounts to raise the questions 
of the gaussianity and statistical isotropy (stationarity) of the 
CMB. New methods of analysis of the CMB must therefore be 
considered. Notably, in addition to the information on scales 
given through a pure spherical harmonics decomposition, the 
scale-space analysis is essential to allow the localization of 
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features on the sky. Such local features might for example 
be associated with non-gaussianity or non-stationarity of the 
statistical distribution from which the CMB arises, or with 
foreground emission such as point sources. First analyses 
on the one-year data of the ongoing Wilkinson Microwave 
Anisotropy Probe (WMAP) satellite mission suggest a 
departure from both statistical isotropy and gaussianity of the 
signal. The various methods applie d exten d from the analysis 
of phase correla tions (|C oles et al. 200 4), N-point correla - 
tion functions (Eriksen et al. 2004; Eriksen et al. 2005), 
bipolar power spectra llHaiian & S ouradeep 2003; 
Haiia net al. 2005t lHaiian & Sou radeep 2 0051). or local 
power spectra (|Hansen et al. 20021 lHansen et al. 2004a ; 
Hansen^£t^L^0j}4bl), to multipole vector s fcopi et al. 2004; 
Katz & Weeks 2 004: Lachiez e-Rev 2004 . In addition, the 
efficiency of the wavelet signal processing for detecting 
non-gaussia nities in the C MB sig nal was also recently es- 
tablished llHobson et al. 1999t iBarreiro & Ho bson 200 ; 
iMartmez-Gonzalez et al. 20021 IStarck et al. 2004 

While no evidence for non-gaussianity was found 
through wavelet analyses in the former COsmic 
Backgroun d Explo rer (COBE) satell ite mission data 
(Barreiro et al. 2000] iCavon et al. 200~fl) . various spher- 
ical wavelet analyses of the one-year WMAP data also 
suggest a departure from non-gauss ianity or statistical 
isotropy dVielva et al. 20041 [ Mukheriee & Wan g 2004t 
McE wen et al. 20051 ICruz et al. 2005^ Other works also 
explicitly established the efficiency of the scale-space wavelet 
processing for point sources detection in the foreground, or 
for technical purp o ses such as denoisin g and deconvo lution 
JSanzetal. 19991: iTenorio et al. 19991: IVielva et al. 2001 
IVielva et al. 2003tlMaisinger et al. 2004 . But a huge amount 



2 



of work is still needed in this context, notably for the identifi- 
cation of not only the position, but also the precise direction, 
and possibly the morphology of the observed fe atures, and 
for their interpretation (see (McEwenet al. 2005) for a first 
approach to the directional wavelet analysis). 

In this context, a new approach to the formalism of spher- 
ical wavelets is developed here. In § [2] we briefly review 
the formalism for the construction of Euclidean wavelets. In 
§ [3J we reintroduce spherical wavelets independently of the 
original group theoretic formalism, in a new, equivalent, and 
self-consistent approach. We adopt a practical philosophy, 
considering wavelets as localized filters which enable scale- 
space analysis, and offer an explicit reconstruction formula 
for the signal considered, from its wavelet coefficients. We 
also prove that the fundamental operation of dilation on the 
sphere is uniquely determined from the requirement of ba- 
sic natural properties. This establishes the uniqueness of the 
wavelet formalism considered. Related technical proofs are 
postponed to appendix [X] In §0] we prove that the inverse 
stereographic projection of a wavelet on the plane gives a 
wavelet on the sphere. It is also established that the stere- 
ographic projection is the unique projection through which 
this correspondence principle holds. This new principle sim- 
plifies the construction of wavelets on the sphere and allows 
to transfer wavelet properties from the plane onto the sphere. 
Notice that it has been mis leadingly suggested in the original 
approach ( Antoine et al. 2002), while not established. The re- 
lated technical proofs are detailed in appendix[2 In § |5] the 
concepts of filter directionality and steerability on the sphere 
are defined from the corresponding notions on the plane. In 
the context of the wavelet formalism, the correspondence 
principle enables to transfer these properties from wavelets 
on the plane to wavelets on the sphere. The property of filter 
steerability allows the computation of the rotation of a filter 
on itself in any direction from a simple finite linear combina- 
tion of basis filters. We also study the angular band limitation 
of steerable filters and explicitly treat the examples of steer- 
able wavelets on the sphere defined as inverse stereographic 
projection of the derivatives of radial functions on the plane. 
In §[5] we discuss the interest of directional and steerable fil- 
ters for the identification and interpretation of local directional 
features in the context of the CMB analysis. A numerical ex- 
ample illustrates our discussion. In §0 we briefly conclude. 

2. WAVELETS ON THE PLANE 

In this section we briefly sketch the well-known formalism 
of wavelets on the plane. 

On the plane as well as on the line, the notion of wavelet 
transfor m of a signal is a powerful me t hod of signal decom - 
position ilTorresani 19951 iMallat 199& lAntoine et al. 2004 . 
We consider here a practical approach for the definition of 
wavelets on the plane, which will easily be translated on the 
sphere. A "mother wavelet" ip(x) is first defined as a localized 
function on the plane, on which affine transformations may 
be applied: translations, rotations, and dilations. Second, the 
wavelet transform of a signal on the plane is defined as the 
correlation of the signal with the dilated and rotated versions 
of the mother wavelet, leading to wavelet coefficients. This 
explicitly defines the scale-space nature of the decomposition. 
In this context, an admissibility condition is finally imposed 
on the mother wavelet by explicitly requiring the exact recon- 
struction formula of the signal from its wavelet coefficients. 
Notice for completeness that, in terms of the original group 
theoretic approach, the wavelet decomposition is defined by 



the construction of the coherent states of the group of affine 
transformations (translations, rotations and dilations) on the 
plane. In this context, a wavelet must satisfy an admissibility 
condition which ensures the square-integrability of the unitary 
and irreducible representations of that group on the Hilbert 
space of square-integrable functions in which the signals are 
defined. This square-integrability implies that the family of 
wavelets obtained by affine transformations from a mother 
wavelet constitutes an over-complete frame in the considered 
Hilbert space, and that an exact reconstruction formula of a 
signal in terms of its wavelet coefficients may be obtained. 
Our more practical considerations lead identically to the same 
wavelet formalism. 

First, the affine transformations are defined as follows. Let 
us consider the Hilbert space of square-integrable functions 
on the plane: g(x) in L 2 (M. 2 ,d 2 x). In the coordinate system 
(o,ox,oy), the point i*on the plane is given in Cartesian co- 
ordinates as x = (x,y), and in polar coordinates as x = (r,(p). 
The invariant measure on the plane, relative to the canonical 
Euclidean metric in R 2 is simply d 2 x = dxdy. Generically, the 
action of an operator on a function in L 2 (R 2 ,d 2 x) is defined 
by the action of the inverse of the corresponding operator on 
R 2 , applied to the function's argument. The operator t (xq) in 
L 2 (R ,d 2 x) for a translation by an amplitude xq = (^o^o) is 
defined in terms of the inverse of the translation t? on points 
in R 2 . Its action reads: 

[t(x Q )g](x) = g(t^x) 1 (1) 

with f? (x,y) = (x + XQ,y + yo). The rotation operator around 
the origin of coordinates r(\) in L 2 (M. 2 ,d 2 x), rotation of the 
wavelet around itself by an angle \ S [0,27r[, is also given in 
terms of the inverse of the rotation r x on points in M 2 . It reads 

[r(x)g](x) = g(r- x 1 x), (2) 

where r x (x,y) follows from the action of the two-dimensional 
rotation matrix r x on the Cartesian coordinates (x,y), or 
equivalently in polar coordinates r x (r,<p) = (r,(p + x)- The di- 
lation operator d (a) on functions in L 2 (R 2 , d 2 x) with a dilation 
factor a G R* is again defined, in terms of the inverse of the 
corresponding dilation d a on points in R 2 . It reads 

[d(a)g](x) = a- 1 g(d- a 1 x), (3) 

with d a (x,y) = (ax, ay) in Cartesian coordinates, or equiva- 
lently d a (r,cp) = (r a (r),ip) for r a (r) = ar in polar coordinates. 
The dilation operator on the plane is uniquely defined (see ap- 
pendix© by the requirement of the following properties. The 
operator d a on the points on R 2 must be a radial (i.e. only af- 
fecting the radial variable r independently of (p, and leaving ip 
invariant) and conformal (i.e. preserving the measure of an- 
gles in the tangent plane at each point of R 2 ) diffeomorphism 
(i.e. a continuously differentiable bijection on R 2 ). The nor- 
malization factor a~ l in Q is also uniquely determined by the 
requirement that the dilation d(a) of functions in L 2 (M. 2 ,d 2 x) 
be a unitary operator (i.e. preserving the scalar product in 
L 2 (R. 2 ,d 2 x), and specifically the norm of functions). 

Second, the analysis of signals goes as follows. The wavelet 
transform of a signal f(x) with the wavelet ip(x), localized 
analysis function in L 2 (M. 2 ,d 2 x), is defined as the correlation 
between the signal f(x) and the dilated and rotated wavelet 
V'x.fl = r (x)d(a)ip, that is as the following scalar product: 

W£(x ,x,a)= / d 2 xipt (x) f (x) 
Jr 2 

= ^Wl/>, (4) 
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for ip? 0lXt a = t(xo)ip x , a . The wavelet coefficients W^(xo,Xi a ) 
represent the characteristics of the signal for each analysis 
scale a, direction x, and position xq. 

Finally, the synthesis of a signal f(x) from its wavelet coef- 
ficients reads: 




W$(xo,x,a) [t(xQ)r(x)L^ipa\ (x). (5) 

In this relation, the operator in L 2 {M. 2 ,d 2 x) is defined by 
the simple division by a constant: [L^g](x) = g(x)/C^. This 
exact reconstruction formula holds if and only if the wavelet 
ip(x) satisfies the following admissibility condition: 

0<C i ,= d 2 k ^\ n <oo, (6) 
Jw \k\ 2 

with the normalization convention ^(k) = J dxe~' k ' x i/j(x) for 
the Fourier transform of functions in the plane. The zero- 
mean is therefore a necessary condition for the wavelet ad- 
missibility in L 2 (M. 2 ,d 2 x): 

J d 2 xi>(x) = Q. (7) 

It is also well-known that under the additional requirement 
that ij)(x) be in L 1 DL 2 (R 2 ,d 2 x), the zero-mean condition (Q 
implies the exact admissibility condition (|6j. Wavelets on the 
plane may therefore be easily built. 

3. WAVELETS ON THE SPHERE 

The formalism of wavelets on the sphere 
was originally establis hed in a group theo- 
retic framework jAntoine & Vanderghevnst 1999t 

Anto ine & Vanderghevnst 19981 ~ lAntoine etal. 2002: 
Demanet & Vanderghevnst 2003; Bogdanova et al. 2005). 
In this section we reintroduce the notion of wavelets on 
the sphere through a new and completely self-consistent 
approach. The resulting formalism is equivalent to the 
original one, but more practical and straightforward. The 
structure of our approach follows, in perfect analogy with the 
formalism of wavelets on the plane introduced in the former 
section. 

For the clarity of the expressions, we denote functions and 
operators on the sphere in uppercase letters, by opposition to 
the lowercase letters denoting functions and operators on the 
plane. Identically to our approach in the plane, a "mother 
wavelet" ^(w) is first defined as a localized function on the 
unit sphere, on which affine transformations may be applied: 
translations, rotations, and dilations. Second, the wavelet 
transform of a signal on the sphere is defined as the corre- 
lation of the signal with the dilated and rotated versions of the 
mother wavelet, leading to wavelet coefficients, defining the 
scale-space nature of the decomposition on the sphere. Third, 
an admissibility condition is imposed on the mother wavelet 
by explicitly requiring the exact reconstruction formula of the 
signal from its wavelet coefficients. We again notice that, 
similarly to the formalism in the plane, in a group theoretic 
approach, the wavelet decomposition is defined by the con- 
struction of generalized coherent states for the affine transfor- 
mations on the sphere (translations, rotations and dilations) 
contained in the conformal group of the sphere, 50(1,3). As 



already stated, our approach leads to the same final wavelet 
formalism. 

First, the affine transformations are defined as follows on 
square-integrable functions G{uS) in L 2 (S 2 ,d£Y) on the unit 
sphere. The point co on the sphere is given in spherical coor- 
dinates as u> = (9,cp). In an orthonormal Cartesian coordinate 
system (o,ox, oy, oz) centered on the unit sphere, the polar an- 
gle, or co-latitude, 9 G [0, 7r] represents the angle between the 
vector identifying uo and the axis oz. The azimuthal, or longi- 
tudinal, angle ip G [0,27r[, not defined though for 8 G {0,7r}, 
represents the angle between the projection of this vector in 
the plane (o,ox,oy) and the axis ox. The invariant measure on 
the sphere, relative to the canonical metric on S 2 induced from 
the Euclidean metric in three dimensions is dVl = d cos 9dip. 
Once again, generically, the action of an operator on a func- 
tion in L 2 {S 2 ,dtt) is defined by the action of the inverse of the 
corresponding operator on S 2 , applied to the function's argu- 
ment. 

The action of a rotation p G 5(9(3) in three dimensions 
on a function G on the sphere, is defined by the operator 
R(p) : G(lo) -> [R(p)G](to) = G(Rr p l Lj) in L 2 (S 2 ,dV). The op- 
erator R p may be decomposed in three consecutive rotations, 
respectively around the axes of coordinates oz, oy, and oz, 
and defined by the Euler angles {<Pq,9q,x), with 9q £ [0,7r] 
and (p ,x G [0,2tt[: R p = R V0 ,e 0tX = R U R k R x- ^ inverse 
rotation R~ l is characterized by opposite Euler angles in the 

reverse order: R^ ofio , x = R -x-e -vo = R-k- x A,k-<po- In the par- 
ticular context of the analysis of functions on the sphere, the 
variable (ipo,9o,x) m the parameter space of the group SO(3) 
may decomposed as (ojq,x) m S 2 <£> [0,27r[, where uiq defines 
a position on the sphere S 2 , and x a direction in [0,27r[ at 
each point. In that context, the rotation operator on a function 
G on the sphere is decomposed as R(p) = R(ojq)R z (x), where 
R(u>o) = R i (tpo)R y (9o) defines the motion or translation of the 
function by ojq, and R z (x) defines its rotation by \ on itself. 
The operator R(ujq) on the sphere for a motion of amplitude 
loq is defined in terms of the inverse of the motion operator 
Ruj on points in S 2 . It reads 

[R(uj )G](u) = G(R-\uj), (8) 

where R^fd^tp) readily follows from the action of the three- 
dimensional rotation matrices R y a and R* on the Cartesian 
coordinates in three dimensions associated with the point oj on 
the sphere. The rotation operator R z (x) of a function around 
itself, by an angle x £ [0,27r[, follows again from the inverse 
of the rotation operator R z x on points in S 2 . It is given as 

[^( X )G]M = g(4" 1 W ), (9) 

where R z x (9, ip) = (9,(p+x) also follows from the action of the 
three-dimensional rotation matrix R* on the Cartesian coordi- 
nates in three dimensions associated with the point u>. 

The dilation operator D(a) on functions in L 2 (S 2 ,dfl), for a 
dilation factor a G M*, is defined in terms of the inverse of the 
corresponding dilation D a on points in S 2 as 

[D(a)G](u) = X 1/2 (a,9)G(D- 1 uj) . (10) 

In spherical coordinates, the dilated point is given by 
D a (0,(p) = (9 a (9),cp) with 9 a (9) = 2arctan(atan((9/2)). This 
corresponds to a linear dilation of tan(#/2): tan{9 a (9)/2) = 
atan(#/2). The dilation operator therefore maps the sphere 
without its South pole on itself: 9 a (9) : 9 G [0, 7r[^ 8 a G [0, 7r[. 
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It is established in appendix [A] that this dilation operator on 
the sphere is uniquely defined by the requirement of the same 
basic and natural properties as on the plane. The dilation D a 
of points on S 2 must be a radial (i.e. only affecting the ra- 
dial variable 9 independently of ip, and leaving cp invariant) 
and conformal (i.e. preserving the measure of angles in the 
tangent plane at each point of S 2 ) diffeomorphism (i.e. a con- 
tinuously differentiable bijection on S 2 ). The conformal factor 
X(a, 9) reads \ l l 2 (a, 9) = a~ l [ 1 +tan 2 (0 /2)] /[ 1 +a~ 2 tan 2 (0 /2)], 
also equal to A'/ 2 (a, 9) = 2a /[(a 2 -l)cosfl+(a 2 +l)l as de fined 
in the original formalism (Antoine & Vanderghevnst 1999). 
The normalization by A 1 / 2 (a, 9) in (II 01 is uniquely deter- 
mined by the requirement that the dilation D(a) of functions 
in L 2 (S 2 ,d£Y) be a unitary operator (i.e. preserving the scalar 
product in L 2 (S 2 ,dil), and specifically the norm of functions). 
Notice that in the limit 9 — > 0, this dilation factor on the sphere 
naturally reduces to first order in 9 to the normalization con- 
stant on the plane: \ l ' 2 (a,6) — ► a -1 . 

Second, the analysis of signals goes as follows. The wavelet 
transform W^{uJo,x,a) of a signal F(co) with the wavelet 
fy(ui), localized analysis function in L 2 (S 2 ,dft) on the sphere, 
is defined as the correlation between F(uf) and the dilated and 
rotated wavelet ^ x . a = R z (x)D(a)^, that is again as the scalar 
product: 



W£(u ,x,a) = J <XlKo,x 



(lo)F(lo) 



with <f 



X-u- 



The 



wavelet 



(11) 
coefficients 



Wy(ido,X) a ) represent the characteristics of the signal for 
each analysis scale a, direction x, and position uiq. 

Third, the synthesis of a signal F(u) from its wavelet coef- 
ficients reads as: 



F(w)-- 



da 



dil 



WU^o,X,a)[R(uj ,x)L^ a ](u;). 



(12) 



In this relation, the operator in L 2 (S 2 ,dST) is defined by 
the following action on the spherical harmonics coefficients 

of functions: L*G /m = G/„,/C^,, with I £ N, m £ Z, and \m\ < 
I. On the contrary to L$ on the plane, the operator does 
not summarize to the division by a constant. Notice that the 
a priori arbitrary choice of scale integration measure c/a/a 3 
is uniquely fixed by requiring the correspondence principle 
discussed in the next section. In appendix|A] a detailed proof 
shows that this exact reconstruction formula holds if and only 
if the spherical harmonics transform ^>i„, of the wavelet fy(ui) 
satisfies the following admissibility condition: 

O^i^IZf" Jl®,,P<». (13) 

for all / G N. This condition is slightly less restrictive 
than the admissibility condition required in the original 
group theoretic formalism for the square-integr abil ity of the 
group representation considered, and given in iB4i . It may 
also be shown that the following relation defines a neces- 
sary condition for the wavelet admissibility in L 2 (S 2 ,dfT) 
(Antoine & Vanderghevnst 1999): 



1+COS0 



(14) 



Finally, in the limit 9^0, the formalism of 
spherical wavelets simply identifies with the formal- 
ism of Euclidean wavelets. This Euclidean limit 
(Antoine & Vanderehevnst 1999) is naturally established 
by direct identification of the relations to @ on the plane, 
one by one with the relations (|8} to J14b on the sphere, to first 
order in 9. In particular, the identification of the admissibility 
conditions (j6jl and Jl 3I > is made obvious in terms of their 
reformulation, as expressed in appendixlHl 

Also notice that the wavelet formalism may be similarly 
defined in the space of integrable functions, on the plane, 
L l (M. 2 ,d 2 x), as on the sphere, /^(S 2 ,^). The dilation opera- 
tor on functions must be redefined consistently in such a way 
that it still preserves the norm of the functions on which it is 
applied. The measure of integration on scales changes and the 
admissibility condition on the sphere is modified consistently. 
In that formalism, the reconstruction formula contains an ad- 
ditional low frequency term, and the corresponding necessary 
admissibility condition on the spher e reduces to an exac t zero- 
mean condition: j s2 dfl^(9,ip) = llAntoine et al. 20021) . But 
we do not develop this formalism here. 

4. CORRESPONDENCE PRINCIPLE 

In the previous sections, we established the wavelet formal- 
ism, independently on the plane and on the sphere. Wavelets 
on the plane are well-known, and may be easily defined in 
terms of the zero-mean condition for a function both inte- 
grable and square-integrable, implying the admissibility con- 
dition (j5Jl. On the contrary, the admissibility condition Jl 3i 
for the definition of wavelets on the sphere is difficult to check 
in practice. In this section, we prove that the inverse stereo- 
graphic projection of a wavelet on the plane leads to a wavelet 
on the sphere. For clarity, the related technical proofs are 
postponed to appendix |B] Beyond its pure theoretical inter- 
est, this new correspondence principle between the wavelet 
formalisms on the plane and on the sphere is of great prac- 
tical use. Indeed, it enables the construction of wavelets on 
the sphere by simple projection of wavelets on the plane. In 
that respect, it also naturally allows to transfer wavelet prop- 
erties from the plane onto the sphere, as illustrated in the next 
section. 

First, appendixlHlestablishes that the projection of a wavelet 
on the plane leads to a wavelet on the sphere under spe- 
cific requirements on the corresponding projection operator. 
That is to say, if the function ip(r,ip) in L 2 (M. 2 ,d 2 x) satisfies 
the wavelet admissibility condition (|6jl on the plane, then the 
function 

*(0,cp)= [irV] d5) 

in L 2 (S 2 ,dft) satisfies the wavelet admissibility condition Jl 31 
on the sphere. The projection operator II between functions in 
L 2 (S 2 ,dfl) and in L 2 (M. 2 ,d 2 x) on the two manifolds is defined 
in terms of the inverse of the corresponding projection oper- 
ator 7r between points on the sphere S 2 and on the plane IR 2 , 
applied to the argument of the function considered. The proof 
of this correspondence principle relies on the requirement that 
the projection operator be a unitary, radial, and conformal dif- 
feomorphism. Let us recall that the unit sphere on which 
spherical wavelets are defined is centered at the origin of 
the orthonormal Cartesian coordinate system (o,ox,oy,oz) in 
three dimensions, with spherical coordinates (9,<p). For sim- 
plicity, and without loss of generality, we consider a geomet- 
rical setting where the plane on which the Euclidean wavelets 
are defined is parallel to the plane (ox,oy), at an arbitrary 
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height zq. The polar coordinates (r, tp) on the plane are de- 
fined relatively to the coordinate system (o,ox,oy). The polar 
angle tp on the plane and the azimuthal angle tp on the sphere 
are therefore identified with one another. We consider specif- 
ically, still without loss of generality, the plane tangent to the 
sphere at the North pole, zo - 1 (see Fig. [Qbelow). With this 
choice of coordinates, the operator tt of projection of points 
from S 2 onto M 2 must naturally be a radial diffeomorphism. 
That is, it must be a continuously differentiable bijection be- 
tween S 2 and R 2 , which only relates the radial variables r on 
the plane and 9 on the sphere independently of tp, and which 
leaves tp invariant. It must also define a conformal mapping 
between the sphere and the plane {i.e. the metric induced on 
the plane by the transformation r(9) is conformally equivalent 
to the Euclidean metric on the plane). This property is essen- 
tial to ensure that the measure of angles and directions (ori- 
entations) is preserved by projection. The projection operator 
II between functions in L 2 (S 2 ,dil) and in L 2 (M. 2 ,d 2 x) on the 
two manifolds must also be unitary (i.e. preserve the scalar 
product between L 2 (S 2 ,dfl) and L 2 (M. 2 ,d 2 x), and specifically 
the norm of functions). 

Second, the correspondence principle may also be extended 
in the following way. The rotation R z (x), acting on functions 
on the sphere through the azimuthal angular variable tp, is 
conjugate to the rotation r(%), acting on functions on the tan- 
gent plane at the North pole through the polar angular variable 
tp, by any radial projection: 



^(x) = n-V( x )n. 



(16) 



As shown in appendix [51 for a unitary, radial and conformal 
projection operator II between L 2 (S 2 ,dfl) and L 2 (M. 2 ,d 2 x), a 
conjugation relation holds between the dilation operators JlQi 
and 0, through the projection II: 



D(a) = U~ 1 d(a)U. 



(17) 



It appears from these two conjugation relations that the oper- 
ations of dilation by a £ R+ and rotation by x £ [0,27r[ of a 
spherical wavelet defined as projection of a Euclidean wavelet 
may be simply performed on the plane before inverse projec- 
tion on the sphere: 

^ x>a (e,tp)= [irVx, J (#,¥>)■ (is) 

The dilated and rotated wavelets on the sphere may therefore 
be built in an extremely straightforward way with the intuitive 
operations of dilation Q and rotation (|2) on the plane, forget- 
ting the corresponding operators (II Oi and on the sphere. 
Also notice that, as all considered operators are unitary, if the 
wavelet on the plane is normalized, the corresponding wavelet 
on the sphere remains naturally normalized. Only motions 
(translations) by luq have to be explicitly performed on the 
sphere as the corresponding operator (|8j is not the conjugate 
of the operator {0 of translation by any xq on the plane. 

Third, appendix |H1 also establishes that the stereographic 
projection is the unique radial conformal diffeomorphism 
mapping the sphere S 2 onto the plane M 2 . The unitary stereo- 
graphic projection operator between functions G in L 2 (S 2 ,dtt) 
and g in L(SH 2 ,d 2 x) and its inverse respectively read 



and 



[nG]« = // 2 (r)G(7r- 1 x), 
[IT 1 *] (u) = tf l l 2 (r(e))g{™). 



(19) 



(20) 



2tan(#/2) and 9(r) = 2arctan(r/2). The diffeomorphism r(9) 
and its inverse 9(r) explicitly define the stereographic projec- 
tion of points on the sphere onto points on the plane, and its in- 
verse. This stereographic projection maps the sphere, without 
its South pole, on the entire plane: r(8) : 9 £ [0,7r[— > [0,oo[. 
Geometrically, it projects a point w = (6, tp) on the unit sphere 
onto a point x = (r, tp) on the tangent plane at the North pole, 
co-linear with uj and the South pole (see Fig. [Q. The con- 
formal factor /x(r) is given as /i I//2 (r) = (1 + (r/2) 2 )" 1 . The 
normalization ^ 1//2 (r) is required to ensure the unitarity of 
the projection operator IT between L 2 (S 2 ,dft) and L 2 (M. 2 , d 2 x). 
Conversely, /i" 1/ ' 2 (r(6')) = 1 +tan 2 (6>/2) is the conformal fac- 
tor associated with the inverse radial conformal mapping ir~ l , 
which ensures the unitarity of IT" 1 . 
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The radial conformal diffeomorphism between points is given 
as %(6,tp) = (r(8),tp) and n~ l (r,tp) = (9(r),tp), and for r(9) = 



FIG. 1. — Relation established between points on the sphere (9,ip) and on 
its tangent plane at the North pole (r, tp) through the stereographic projection 
it and its inverse 7r _1 . The same relation holds between functions living on 
each of the two manifolds, as illustrated by the shadow on the sphere and 
the localized region on the plane. The corresponding operators II and II -1 
enable the establishment of a correspondence principle between the wavelet 
formalisms on the plane and on the sphere. 

Let us consider one example, as illustration of the cor- 
respondence principle through the stereographic projection. 
On the plane, the derivatives of gaussians in a specific di- 
rection, say x, are well-known examples of wavelets. The 
normalized (negative) first derivative of gaussian in direction 
x reads, after dilation by a and rotation by X- 4 , x,a( r t ( fi) = 
yjlj-x re~ rl l 2a ~ cos(tp - x) / a 2 ■ The corresponding wavelet on 
the sphere, normalized, dilated by a and rotated by x> but 
still at the North pole, is simply obtained by the action of 
the inverse stereographic projection (I20> : ^ x a (9,tp) = (1 + 

tan 2 (6>/2)) x /27^2tan(0/2)<T 2tan2(e ' /2) / fl2 cos(tp - x)/a 2 . No- 
tice that, depending on the application, the exact zero mean of 
a filter is often an appreciated property, as the corresponding 
filtering completely erases constant signals. In the L 2 formal- 
ism, this property may be considered as an additional practical 
requirement for the choice of the mother wavelet. As illus- 
trated by the present example, the stereographic projection of 
any odd-numbered derivative of a radial function already sat- 
isfies this zero mean condition, in addition to the admissibility 
condition d 1 3I > . 

Finally, these developments necessitate the following re- 
marks. Notice on the one hand, that the conformal rela- 
tion between the Euclidean invariant measure on the plane 
and the measure induced from the stereographic projec- 
tion, r(9)dr(9)dtp = /i" 1 (r(9)) sin OdOdtp with r(9) = 2 tan(0/2), 
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readily implies the correspondence of the necessary wavelet 
admissibility conditions on the plane and on the sphere. If a 
function satisfies the necessary wavelet admissibility condi- 
tion on the plane, then its inverse stereographic projection 
satisfies the necessary wavelet admissibility condition dl4t on 
the sphere. This result gives us intuition on the otherwise non- 
trivial complete demonstration of the correspondence princi- 
ple ill 5b established in appendix[B] which links the necessary 
and sufficient wavelet admissibility conditions. On the other 
hand, we also emphasize that, in the limit 9 — » 0, we get the 
first order relation r(9) — > 9, which corresponds to the identi- 
fication between the considered portion of the sphere around 
the North pole, with the tangent plane at that same point. The 
factors of unitarity also tend to unity. The stereographic pro- 
jection summarizes to the identity operator. The wavelets on 
the plane are therefore identified with their projection on the 
sphere, in complete coherence with the Euclidean limit dis- 
cussed in §|3] 

5. FILTER DIRECTIONALITY AND STEERABILITY 

In this section, we discuss the notions of directionality and 
steerability of filters on the plane and on the sphere. On the 
sphere, these developments constitute a new advance for the 
scale-space analysis of signals. Their interest will be illus- 
trated in the next section in the context of the CMB analysis. 
In a first subsection, we recall the well-known notions of di- 
rectionality and steerability on the plane. We introduce the 
corresponding definitions on the sphere. We also show that 
these notions are transferred from the plane onto the sphere 
through the stereographic projection. In the context of the 
wavelet formalism, the correspondence principle established 
in the former section therefore enables to transfer the direc- 
tionality and steerability of wavelets, from the plane onto the 
sphere. In a second subsection, we characterize steerable fil- 
ters in terms of their band limitation in the Fourier index con- 
jugate to the angular variable ip. In the last subsection, we fi- 
nally develop the examples of steerable wavelets on the sphere 
defined as inverse stereographic projection of the derivatives 
of radial functions on the plane. 

5.1. Definitions and correspondence 

First, we recall the notions of directionality and steerability 
on the plane. Locally, at each point of coordinates x = (r, ip) 
on the plane K 2 , directions are defined in the tangent plane, 
in terms of the rotation angle x G [0,27r[. The origin of an- 
gles (x = 0) in the tangent plane is defined by the direction 
tangent to the line passing through the point considered, from 
the origin o of coordinates, and making angle ip with the axis 
ox (direction of increasing r) . The directionality of a filter 
g(r,ip) in L 2 {M. 2 ,d 2 x) on the plane may be measured through 
its auto-correlation function, defined as the scalar product of 
the rotations of the filter in two different directions x an d x'> 
and depending on the difference A\ = X~x' : 

CHA X )=(r(x)g\r(x / )g)- (21) 
The steerability of a directional filter g(r, <&) is d efined 
in L 2 (R 2 ,d 2 x) on th e plane jFreeman & Adelson 199 It 
ISimoncelli et al. 19921) by the requirement that any rotated 
version r(x)g of the filter be expressed as a linear combina- 
tion of a finite number of rotations of the filter in specific di- 
rections Xm- This is mathematically defined by the relation 
[r(x)g](r,<p)= Y,m=i k >n(x)MXm)g](r,<p), where the weights 
km(x) > with I < m < M, and M £ N, are called interpola- 
tion functions. More generally, one requires that r(x)g be ex- 
pressed in terms of a finite number of independent basis filters 



gm(r,tp) in L 2 (R 2 ,d 2 x), which are not necessarily specific ro- 
tations of g(r, ip): 

M 

[r(x)g] 0, f) = ^2 k >" ^Sm(r, <p) ■ (22) 

m= 1 

Second, we can define the notions of directionality and 
steerability on the sphere in perfect analogy with their def- 
inition on the plane. Locally, at each point of coordinates 
u> = (9, ip) on the sphere S 2 , directions are defined in terms 
of the third Euler angle x € [0,27r[, which identifies the di- 
rections in the tangent plane. The origin of angles (x = 0) 
in the tangent plane is defined by the direction tangent to the 
meridian passing through the point considered (direction of 
increasing 9). Let us recall that the non-existence of differ- 
entiable vector fields on S 2 rules out the definition of direc- 
tions globally on the sphere. The directionality of a filter 
G(9,ip) in L 2 (S 2 ,dn) on the sphere is measured through its 
auto-correlation function, also defined as the scalar product 
of the rotations of the filter in two different directions x an d 
x', and depending on difference Ax = x~X >: 

C c (A x )=(R i (x)G\R i (x')G). (23) 

The steerability of a directional filter G(9,ip) is defined in 
L 2 (S 2 ,dfl) on the sphere by the relation 

M 

[R* ( X ) G](6,tp)=J2 k m (X) G m (6,<p). (24) 

m=l 

Again, the weights k m (x), with 1 < m < M, and M £ N are 
the interpolation functions. The basis filters G m (9, ip) are not 
necessarily specific rotations of G(9,(p). 

Third, we show that these properties are transferred from 
the plane onto the sphere through the stereographic projec- 
tion. Considering the stereographic projection G = IT'g, the 
unitarity of the operator II ensures that the auto-correlation 
function of the projected filter is identical to the auto- 
correlation function of the original filter on the plane: 

C (n_lg) (Ax) = C«(Ax). Moreover, if the relation of steerabil- 
ity 1221 holds for g(r,ip) on the plane, then the relation of 
steerablity 1241 obviously holds for G = U~ l g on the sphere, 
with the same interpolation functions k m (x) > f° r 1 < m < M, 
and M £ N, and for the basis filters G„,(9,(p) defined as 
G,„ = H~ l g,„. This property explicitly relies on the conjuga- 
tion relation d!6l > induced by any radial projection operator. 
If wavelet filters are considered, the correspondence princi- 
ple established in §|4] therefore enables to transfer the direc- 
tionality and steerability of wavelets, from the plane onto the 
sphere. 

We finally discuss the interest of the concepts introduced, 
on the plane as on the sphere. Let us consider on the one 
hand the notion of filter directionality. Any non-axisymmetric 
filter will be defined as directional in the sense that its auto- 
correlation is not a constant function of Ax- We consider as a 
good directional filter, a filter for which the auto-correlation, 
J2 1 i or d23l >. is a rapidly decreasing function of Ax- In 
this regard, the ideal directional filter would have the ex- 
pression of a S distribution in the angle ip, in such a way 
that its auto-correlation be a 6(A\) distribution. Notice that 
other definitions of directionality on t he plane and on the 
sphere may be found in the literature jAnto ine et al . 2002t 
Demanet & Vanderghevnst 2003; Antoine et al. 2004). Our 
definitions have the non-negligible advantage of correspond- 
ing to one another through the stereographic projection be- 
tween the sphere and the plane. In these terms, the peakedness 
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of the auto-correlation function is a measure of the sensitivity 
of the filter to directions. From a practical point of view, it 
also measures how the identification of directions, in terms of 
the maximization of the filtering coefficients of a considered 
signal, is sensitive to any kind of noise inevitably affecting the 
estimation of these coefficients. The directionality is conse- 
quently a key criterion for the choice of the filter for the iden- 
tification of local directions of a signal. Consider on the other 
hand the notion of filter steerability. Through a linear filter- 
ing, if a relation of steerability holds for a filter, i22\ or (I24> . 
then the same relation holds for the filtering coefficients of 
the signal considered. This is specifically true in the analysis 
of a signal with wavelet filters. The steerability therefore al- 
lows the computation of filtering coefficients of a signal in all 
directions at the cost of the computation of M filtering coef- 
ficients. This property may therefore reduce the computation 
cost of filtering by a non-negligible factor. Its direct interest 
in the perspective of the CMB analysis is suggested in § [6] 
In the following subsection, we discuss the band limitation in 
the Fourier index conjugate to the variable tp for directional 
and steerable filters. 

5.2. Angular band limitation 

We show here that the notions of ideal directionality and 
steerability represent competing concepts on the plane as on 
the sphere, in terms of the band limitation of the considered 
filters in the Fourier index conjugate to the angular variable tp. 
The following results are discussed on the plane, for functions 
g(r,tp) in L 2 (R. 2 ,d 2 x). They may be identically read on the 
sphere for G{6,tp) in L 2 (S 2 ,d£Y), through the substitution of g 
by G. Let us thus consider the Fourier decomposition of the 
filter g(r,tp) in the variable tp: g(r,<p) = Y^H-NSni^e"^ /2n. 
The function g„(r) stands for the n' h Fourier coefficient, and 
N + \ represents the band limitation in the Fourier index n: 
n <N+\. On the one hand, an ideal directional filter has an 
angular dependence associated with a 5(tp) distribution, which 
corresponds to a null angular width, that is no band limit, 
N — > oo, and constant Fourier coefficients. On the other hand, 
conceptually, the basis filters of a steerable filter must have a 
non-zero angular width (hence also the filter itself by linearity 
of the steerability relation). This ensures that they are sensi- 
tive to a whole range of directions. In that case only, one may 
imagine to steer the filter in all directions from a finite num- 
ber, M in relation J22t . of these basis filters. The non-zero 
angular width is naturally associated with a band limitation N 
of the steerable filter and its basis filters. It may be proved rig- 
orously that, if T is the finite number of non-zero coefficients 
g„(r) for a filter g(r, tp), then this filter is steerable and T is the 
minimum number of basis filters g m (r, tp) required in relation 
d22t to steer g( r, tp). In other words, the fo llowing inequality 
holds: M >T dFreeman & Adelson 19911) . Consequently, if 
the filter g(r, tp) is steerable with a number M of basis filters, 
then it has a finite number T of non-zero coefficients and, as 
suggested here above, it is inevitably limited in band at some 
band limit N. The notion of ideal steerable filter (M small, 
hence N finite) is therefore clearly in opposition with the no- 
tion of ideal directional filter (N — > oo, hence M — > oo). 

5.3. Examples of steerable wavelets 

The present subsection introduces the inverse stereographic 
projection of derivatives of radial functions on the plane as 
examples of steerable wavelets on the sphere. On the plane, 
if a sufficiently regular radial function cj>{r) is considered, its 



N th derivative in direction x, 

^(r,^) = ^[0(r)], (25) 
satisfies the wavelet admissibility condition for any N > 1 and 
may be called a wavelet. In terms of the Fourier index n conju- 
gate to the angular variable ip, this function is limited in band 
at the level N + 1 . It contains T = N + 1 non-zero coefficients 
of indices n = -N, -(N-2), 0, N-2, N for N even, and 
n = -N, -(N-2), 1,-1, ...,N-2,NfovNodd. The rotation 
by an angle \ of the N th derivative in direction x of a radial 
function is given as the N' h derivative in direction u(x) = r x x 

It may thus be expanded as: 



:#(^) = ^ x) [</>(r)]. 



2 

E 

l •■■'«=! 



Mi,...^/'''!'"^ (r,tp) 



(26) 



with x\ = x and xj = y, and where the coordinates of u(x) 
read (1*1,1*2) = (cosx,sinx). It is therefore a steerable fil- 
ter, expressed in terms of M = 7^ = N +1 = T combinations 

with repetitions of the basis filters tp di °> = df~ k d~[<M.r)], for 
0<k<N. In the following paragraphs, we explicitly consider 
the examples of the first and second derivatives of a gaussian. 

The first derivative of a radial function in direction i on the 
plane reads: 

tp 9 ' (r, >p) = d r (f>(r)cos tp, (27) 

where d r stands for the radial derivative. It therefore has an 
angular band limit at + 1 = 2. It is far from being an ideal 
directional filter as its auto-correlation function reads for a 
normalized filter: C 9x (Ax) = cos(Ax). In terms of steerability 
however, a first derivative is an ideal filter, as it only requires 
N+l = 2 weights. The steerability relation reads, in terms 
of the specific rotations ip 9s and tp 9 ' at x = and x = 7r /2 
respectively: 

</, tp) = tp 9s (r, tp) cos x + ip 9i (r, tp) sin x- (28) 
As a concrete example, let us consider once more the normal- 
ized (negative) first derivative of a gaussian already discussed 
at the end of §|4] It is given here through the general relation 
J27> for <j)(r) = -^2/ite~ r I 2 . The section of this wavelet at 
y = is illustrated by Fig. while its section at constant* has 
a gaussian shape: 



FIG. 2. — Section at y = of the first derivative of gaussian in direction x 
on the plane. 



The inverse stereographic projection on the sphere of the 
wavelet rotated by x ar, d dilated by a (see relation (I18» . is 
given, as already discussed, by: 



d s (gauss) t 



2 2 



(9,tp) = J-^ 1+tan 2 - 



tan-cos(^-x) 



-2tan 2 (S/2)/a 2 



(29) 



8 



Fig. [3]illustrates the directionality of first derivatives of radial 
functions and their steerability relation d28l >. understood on 
the sphere, through the specific example of the first gaussian 
derivative: 






FIG. 3. — First derivative of gaussian on the sphere for a dilation fac- 
tor a = 0.4: (top left) flj, (top right) dy, and (bottom) dc,( x ) with \ = t/4. 
Dark and light regions respectively identify negative and positive values of 
the functions. 



The second derivative of a radial function in direction x on 
the plane reads: 



o2 \ O O 

ip * (r, ip) = -d r <j>(r)sm cp+d,(f>(r)cos <p. 
r 



(30) 



It is limited in band at + 1 = 3. Its auto-correlation func- 
tion reads C a ? ( Ax) = A + B cos 2 A\, with the values A and B 
functions of (f>(r), and A + B = 1 for a normalized filter. A sec- 
ond order derivative is therefore not necessarily a better direc- 
tional filter than a first order derivative as its auto-correlation 
function is not generically better peaked. The general steer- 
ability relation (122> holds, the rotated filter being expressed in 
terms of basis filters, which are not specific rotations of the 
second derivative itself. It requires N+ 1 = 3 weights: 



iht (r,<P) = (r, f) cos 1 x + i> d * (r, <fi) ^ X 
+il) a * d *(r,ip) sin2x- 



(31) 



Again, we consider the example of the normalized (neg- 
ative) second derivative of a gaussian, that is for </>(r) = 

-^4/3vre-'"/ 2 in relation (I30i . The section of this wavelet 
at y = is illustrated by Fig. |4] while its section at constant x 
has a gaussian shape: 





A 


^a|(gei»»») 







FIG. 4. — Section at y = of the second derivative of gaussian in direction 
x on the plane. 



The inverse stereographic projection on the sphere of the 



wavelet rotated by x an d dilated by a reads 
a V jtt \ z 



2tan 2 (9/2)/a 2 



1 — j tan - cos (<p—x) 



(32) 



Fig.|5]illustrates the directionality of second derivatives of ra- 
dial functions and their steerability relation ( 13 1> . understood 
on the sphere, through the specific example of the second 
gaussian derivative: 

z z z 

x ^^^^^ y x ^^^^^ y x ^^^^^ y 

z 

^^^^^ 



FIG. 5. — Second derivative of gaussian on the sphere for a dilation factor 
a = 0.4: (top left) c9?, (top centre) c9?, (top right) dgdf, and (bottom) 8^„j 
with x = 7r/4. Dark and light regions again respectively identify negative and 
positive values of the functions. 

6. CMB LOCAL DIRECTIONAL FEATURES 

In this last section, we discuss the interest of directional 
and steerable filters for the identification and interpretation of 
local directional features on the sphere, with direct application 
to the CMB analysis. First, we show that steerable filters may 
efficiently detect local directional features on the sphere with 
the same angular precision as ideal directional filters. Second, 
the discussion is illustrated by a simple numerical example. 
Third, we emphasize on how the property of filter steerability 
is essential to reduce the computation cost in the search for 
directional features through the wavelet analysis of the CMB. 

The theoretical angular resolution power of any directional 
filter remains infinite, independently of its possible steerabil- 
ity. This infinite resolution assumes however that the mor- 
phology of the signal considered is known a priori, in partic- 
ular in the variable <p. In the context of the CMB analysis, 
this might be the case when looking for the imprint of cosmic 
strings or other pre-defined features in the background radi- 
ation. In that case indeed, the wavelet coefficient is known 
analytically as a pure function of the difference A\ = X~X* 
between the wavelet rotated by x< an d tne direction x* of 
the considered feature. This direction x* ma Y therefore be 
identified exactly in the analysis of the wavelet coefficients 
of the signal. Notice however that for a given experimen- 
tal precision, the resolution power of filters is directly related 
to their directionality, through the peakedness of their auto- 
correlation function, and is therefore a function of the angular 
band limitation and the number of weights M in the case of 
steerable filters. 

The following simple example illustrates the theoretical in- 
finite angular resolution power of steerable filters. 
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FIG. 6. — Test-signal: elongated feature centered at the point u>* = 
(6*,tp*) = (7r/2,7r/2) of the sphere, and making an angle \* = n/3 with 
the meridian at that point. 
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FIG. 7. — Analysis-wavelet: second gaussian derivative in direction x, 
projected on the sphere. It is shown centered at the point lu* = (9* ,<p*) = 
(7r/2,7r/2) of the sphere, at a scale a = 0.05, and making an angle Xo = X* + 
7r/2 = 5-7I-/6 with the meridian at that point, that is in the configuration which 
maximizes the corresponding wavelet coefficient of the test-signal. 

First, the test-signal is defined as an elongated feature 
centered at the point u>* = (0*,cp*) = (7r/2,7r/2) of the 
sphere, and making an angle \* = n /3 with the merid- 
ian at that point. The signal is analytically defined by the 
function F(9,(p) = exp[-(cosx*jc-sinx*z) 2 /2cr 2 ]exp[-16(3'- 
l) 2 ], with the Cartesian coordinates of points of the 
sphere related to their spherical coordinates by (x,y,z) = 
(sin 9 cos Lp, sin 9 sin <p, cos 9), and for a half width a = 0.05 
(see Fig. Second, the signal is analyzed by its correlation 
with a second derivative of a gaussian iff 9 ', that is in terms of 
its wavelet coefficients Wg 2 (coo,X^ a ) at eacn point u>o of the 

sphere, and for any scale a and direction \ of the wavelet. 
We might have equivalently chosen the first gaussian deriva- 
tive or any other steerable wavelet. Notice that the zeros of 
the second gaussian derivative at the North pole are located 
at f?o = 2 arctan(a /2) in the direction of the wavelet ip = y, 
and in the opposite direction ip = x + 7T ( see equation ( 1321 . 
and Fig. yj). This angular opening 6q may be understood 
as a qualitative measure of the half width of the wavelet. For 
a dilation factor a = 1, we get 9q ~ 0.9rad. For a small di- 
lation factor a, the half width is simply given as 9o ~ arad. 
We consider a highly localized analysis-wavelet for a dila- 
tion factor a = 0.05, which corresponds to a typical half width 
9q ~ 0.05 rad. It is therefore of the same size as the width of 
the feature in the direction \* + -k j2 as a = a = 0.05, but much 
smaller than its elongation in the direction \* ■ At the North 



pole, the non-rotated second gaussian derivative with a = 0.05 
has essentially zero mean in direction i (this is only exact in 
the Euclidean limit, that is a — > 0), and a maximum mean in 
direction y. At any point u>q on the signal, the wavelet trans- 
form is therefore minimum if the wavelet is directed along the 
signal, which is essentially constant on an interval of the size 
of the wavelet. On the contrary, it is maximum if it is directed 
at an angle ir/2 relative to the feature, direction in which it 
has the same scale as the signal itself (see Fig. [7}. 

Through the correspondence principle and the linearity of 
the wavelet filtering, the steerability relation i31\ for the sec- 
ond derivative of a gaussian may be equivalently written for 
the spherical wavelet coefficients as: 

W # (X) = W k COs2 X + W (? Si " 2 X + W V t Sin 2x ' (33) 

where we dropped the dependence of each coefficient in loq 
and a. Each basis coefficient is evaluated at x = for the 
corresponding basis wavelet. This relation may also be writ- 
ten in a form similar to the expression of the auto-correlation 
function of the second gaussian derivative. This is natural 
as the auto-correlation function is nothing else but the par- 
ticular wavelet coefficient of the rotated wavelet analyzed 
by itself. We get indeed Wg 2 (x) = ^'+#'cos2(x-Xo), for 

A' = and B' = (W^ -W£)/2cosxo, with a max- 

imum at xo defined by 

tan2Yo= — s — —e- (34) 

Or ai 

At each point ujq considered, the cost of the analysis is there- 
fore reduced to the computation of the three wavelet coef- 
ficients for the basis wavelets at the chosen scale a = 0.05, 
and in their original direction \ = 0. The direction x* of the 
feature at that point is given as \* = Xo~ 7T /2- Fig- 13 rep- 
resents the wavelet coefficient W F ^ (x) as a function of x, re- 
s- 
suiting from this directional analysis at the point uj$ = ui* = 
(7r/2,7r/2), through the relation fl33t . 
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FIG. 8. — Wavelet coefficient as a function of \ for the considered 
elongated feature, analyzed with a second gaussian derivative at the point 
w o = (t/2, 7t/2) and at the chosen scale a = 0.05. 

The direction of the feature is precisely recovered at x* = 
7r/3, up to negligible numerical errors. This clearly illustrates 
the fact that steerable wavelets have an infinite angular resolu- 
tion power for an ideal experiment where no error affects the 
signal. 

In conclusion, the steerability of wavelets is a property of 
fundamental interest in the analysis of local directional fea- 
tures on the sphere. First, the filter steerability allows a 
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drastic reduction in computation cost for analyzing all pos- 
sible directions of features at each point on the sphere, as 
only a small number M of basis directions must be con- 
sidered at each point. Second, as the above example illus- 
trates, this reduction is achieved theoretically without any 
precision loss in the identification of the directions. In the 
context of the CMB study, let us recall that local direc- 
tional features may be associated with non-gaussianity, sta- 
tistical anisotropy, or foreground emission. Their identifica- 
tion requires for each analysis scale a, the correlation of a 
filter in all directions \ an d at a U points uq of the sphere 
with numerous theoretical simulations of the signal. Such 
an analysis is currently hardly affordable in terms of com- 
putation time at the already high resolution level of the CMB 
maps, when one wants the same precision in the identification 
of the direction (> 10 3 sampling points) as of the localiza- 
tion (10 6 sampling points on the sphere) of singular features 
(McEwen et al. 2004; McEwe n et al. 20051) . The property of 
filter steerability may therefore be used to reduce the com- 
plexity of such calculations and eventually render them acces- 
sible. However, Fig. |8]also illustrates the non optimal peaked- 
ness of the wavelet coefficients of the considered signal as a 
function of directions (W£(x) = A' + B' cos2(x~ xo)), equiv- 

alently to the peakedness of the wavelet auto-correlation. As 
already discussed, in a practical situation, a better directional- 
ity of the filter would enhance the stability of the measurement 
of directions relatively to noise. The specific choice of filter, 
optimizing the compromise between ideal directionality and 
steerability, will depend very much on the application consid- 
ered, in terms of noise, precision requirements, and computa- 
tional resources. 

7. CONCLUSION 

The recent developments in the analysis of the cosmic mi- 
crowave background (CMB) radiation ask for new methods 



of scale-space signal analysis on the sphere, notably for the 
identification of foreground emission, or for testing the hy- 
potheses of gaussianity and statistical isotropy of the CMB, 
on which the concordance cosmological model relies today. 

In the work presented here, we reintroduced the formalism 
of wavelets on the sphere in a practical and self-consistent 
approach, simply understanding wavelets as localized filters 
which enable a scale-space analysis and provide an explicit 
reconstruction of the signal considered from its wavelet co- 
efficients. We also proved a correspondence principle which 
states that the inverse stereographic projection of a wavelet on 
the plane gives a wavelet on the sphere, and allows to transfer 
wavelet properties from the plane onto the sphere. In that con- 
text, we finally defined and discussed the notions of direction- 
ality and steerability of filters on the sphere for the analysis of 
local directional features in the signal considered. 

The practical formalism introduced provides a method for 
the detection of local features in the CMB, and the iden- 
tification of their precise direction, through the analysis of 
the wavelet coefficients of the signal. Notice however that 
these generic developments of signal processing on the sphere 
may find numerous applications beyond cosmology and astro- 
physics, as soon as the data to be analyzed are distributed on 
the sphere. 
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APPENDIX 

TECHNICAL PROOFS FOR WAVELETS ON THE SPHERE 

In this first appendix, we establish technical proofs for the wavelet formalism on the sphere defined in §13 First, we prove the 
uniqueness of the unitary, radial, and conformal operator D(a) of dilation on functions in L 2 (S 2 ,d£Y) on the sphere, and give its 
expression. Second, we explicitly prove the admissibility condition for wavelets on the sphere (1131 . This condition is required to 
ensure the exact reconstruction formula (I12> for any signal F(ui) from its wavelet coefficients il It . in a decomposition with the 
wavelet in L 2 (S 2 ,dft). 

Uniqueness of the dilation operator 

The dilation operator D(a) on functions in L 2 (S 2 ,dil) is defined in terms of the inverse of the corresponding dilation D a on 
points in S 2 , applied to the argument of the function considered. The dilation operator D a must be a radial and conformal 
diffeomorphism. A radial operator only affects the radial variable 9 independently of <p, and leaves <p invariant. This property 
is natural to define a dilation of points on the sphere relative to the North pole. In this context, the operator D(a) on functions 
in L 2 (S 2 ,dft) takes a generic form given by relation (llOi . in terms of a radial diffeomorphism D a (9,(f) = {8 a {8),(p) on S 2 . The 
diffeomorphism a (0) and the function A'/ 2 (a, 9) have to be determined from the additional properties required. 

First, the diffeomorphism D a may be understood as a coordinate transformation on the sphere, D a : u> = (x l ,x 2 ) = (8,<p) — > 
u' = (x n ,x' 2 ) = (8 a ,(p) G S 2 . It must be conformal in order to preserve the measure of angles and directions, which are 
defined locally in the tangent plane at each point of S 2 . This explicitly requires that the transformed metric g' t j(9 a ,(p) = 

(dx k I 'dx")(dx' j dx'-')gki(9(6 a ),ip) is conformally equivalent to the original metric gy(0 a ,</j) on the sphere, with i,j,k,l G {1,2}. 
We consider the metric induced on the sphere from the Euclidean metric in three dimensions: gij(9, ip) = diag(l , sin 2 9). The con- 
formal equivalence reads by definition: g^(0 a ,<p) = e^ < - a ' e ^ e "^gij(9 a ,ip), for some strictly positive conformal factor e^ {a - e \ that is 
for a real function <fi{a,9). This condition implies that the operator D a is linear in tan 9/2: X.w9 a {9)/2 = a(a)tan9/2. The function 
a(a) must be strictly increasing in the dilation factor a G R+, with the limits a(Q) = and a(oo) = oo. The group structure for the 
composition of dilations must also be fulfilled. This constrains a(a) to a(a) = a a ° for a strictly positive exponent ao G R+. We 
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take, without loss of generality, a function linear in the dilation factor a, corresponding to cxq = 1 : 

tan = atan-. (Al) 

2 2 

Any other choice for oq would simply correspond to a rescaling of the dilation factor. By dilation, the sphere without its South 
pole is thus mapped on itself: 9 a {8) : 9 6 [0,7r[— » 9 a G [0,7r[. The conformal equivalence also determines the conformal factor: 
e -4>(afi)li = fl[1 +t an 2 (6l/2)]/[(l +a 2 tan 2 (6»/2))]. 

Second, the unitarity of the dilation operator D(a) on functions in L 2 (S 2 ,dil) identifies X(a,9) with the conformal factor: 
X(a,0) = e^v^ = e-M" 1 ' 6 ), or equivalently A 1 / 2 (a,6>) = a' 1 [I +tan 2 (6»/2)]/[(l + fl- 2 tan 2 (6»/2))]. This identity simply rehes 
on the relation between the invariant measure on the sphere and the measure obtained after the conformal transformation D a : 
sm9 a (9)d9 a (9)d(p = e-^ a > 0> sm9d9dip. 

On the plane, the construction of the dilation d(a) on functions, unique unitary, radial, and conformal operator in L 2 (M. 2 ,d 2 x), 
relies on an identical reasoning. The canonical Euclidean metric on the plane in polar coordinates simply reads gij(9,(p) = 
diag(l,r 2 ), for i,j 6 {r,ip}. The expression (0 trivially follows. 

Establishment of the admissibility condition 

First, we recall the expressions for the Fourier decomposition of a function G(uS) in L 2 (S 2 ,d£Y) on the sphere. The spherical 
harmonics F/ m (u>), with / G N, m G Z, and \m\ < I, form an orthonormal basis in L 2 (S 2 ,d£l): 

dQ Y,* m (u) Y Vm , (u) = 6„, <W . ( A2) 

s 2 

The direct and inverse spherical harmonics transforms of a function G(u>) are respectively defined as: 

G lm = [ dnYf m (u)G(u) (A3) 



JS 2 

G{uo) = Y J Y.^' mYlm(uo) - (A4) 

im \m\<l 

Second, we introduce the transformation law of functions on the sphere under rotation, in terms of the Wigner D-functions 
jBrink & Satchler 19931) . Let p be an element of 50(3), with p = (ip,9,x), in a decomposition in the Euler angles tp, 9, and 
X- The Wigner D-functions D l nm (p), with / 6 N, m,n 6 Z, and \m\, \n\ < I, are the matrix elements of the irreducible unitary 
representations of weight I of the rotation group in the space of square-integrable functions L 2 (SO(3),dp) on S0(3), with the 
invariant measure dp = dipdco&9dx- By the Peter- Weyl theorem on compact groups, the matrix elements D l * n also form an 
orthogonal basis in L 2 (SO{3),dp), with the following orthogonality relation: 

/ dpD' mn (p)D l J n ,(p)= — — -<5//<<W<W. (A5) 

Let us consider the decomposition p = (too, x), with uiq = (6o,<po) identifying a point on the sphere S 2 , and \ £ [0, 27r[ identifying 
a direction at each point. The action of the corresponding operator R(u>o,x) on a function G(uf) in L 2 (S 2 ,d£l) on the sphere reads 
in terms of its spherical harmonics coefficients and the Wigner D-functions: 

[R U^X) G\ m = D' m „ (wo , X ) Gi n ■ (A6) 

\n\<l 

Finally, we prove the admissibility condition for a wavelet on the sphere. Let F(u>) be a function in L 2 (S 2 ,dil) on the sphere, 
and fy(oj) in L 2 (S 2 ,d£Y), the wavelet considered for the decomposition into the coefficients Wy(u>o,x,a) give n in il 1> . We w ant 
to establish the condition under which the explicit reconstruction formula (I12> holds. From the expressions (IA4> and (IA6> and 
the definition of the operator L*, the function [R(uj(),x)Ly , i l a ](uj) = [L$^a](R^} x <^>) takes the form 

[R(LJ ,X)L^aKu) = Y / Yl ^ D L(UO,X)($a)lnYlm(u)- (A7) 
/eN|m|,|»|</ * 

The wavelet coefficient Wy(uJo,Xi a ) defined in (lilt may be written as: 

w£(LJo,X,a) = J2 £ D l :„(u;o,x)(^)lF, m . (AS) 
ten M,|n|<* 

Inserting these last two expressions in dl 2i . and using the orthogonality relation ( IA5t for the Wigner D-functions, we obtain: 



F(u) = YY,F lm Y lm ((j)-±- 

l£N\m\<l * 



87r 2 f + °° da , 2 



21+1 ^ 

\n\<l 



(A9) 



From this last expression it is obvious that the reconstruction formula H2\ holds if and only if the coefficients defined in Jl 3i 
are finite and non-zero for any / G N. This explicitly establishes the wavelet admissibility condition Jl 3I > on the sphere. 



12 



TECHNICAL PROOFS FOR THE CORRESPONDENCE PRINCIPLE 

In this second appendix, we prove the correspondence principle defined in § |4] which states that the inverse stereographic 
projection of a wavelet ip(x) in L 2 (M. 2 ,d 2 x) on the plane, thus satisfying the admissibility condition (jfji, gives a wavelet ^(w) in 
L 2 (S 2 ,dQ) on the sphere, satisfying the admissibility condition Jl 3i . First, we reformulate the admissibility conditions on the 
plane and on the sphere. Second, we establish explicitly that the projection of a wavelet on the plane gives a wavelet on the sphere, 
in terms of these reformulated admissibility conditions. In that respect, we only require that the projection n between functions 
in L 2 (S 2 ,dn) and in L 2 (M. 2 ,d 2 x) be a unitary, radial, and conformal diffeomorphism. Finally, we prove that the stereographic 
projection is the unique projection operator satisfying these properties. 

Reformulation of the admissibility conditions 

First, let ip(x) be a function in L 2 (R 2 ,d 2 x). We show that the wavelet admissibility condition l|6} on the plane applied to ip(x) is 
equivalent to the condition: 

pin p+co ^ a p 

0<ll= d X — / d 2 xoK^, x , a |/>| 2 <oo, (Bl) 
Jo Jo a Jr 2 

for any f(x) ^ in L 2 (R 2 , d 2 x). To prove this statement, we simply notice that the scalar product is preserved up to a factor 2ir by 
the Fourier transform. This implies that 27r(-0 i - o , x , a \f) = ((4>x ,x-a)\f)> w i m (i>x ,x,a)(k) = aijj{af^k)e~ lk ' Xa . Using the orthogonality 
relation of the imaginary exponentials, the integral therefore takes the form 



p2tt f+oo , 

Jo Jo a 



(B2) 



where 1 1/| | stands for the norm of f(x) in L 2 (R 2 , d 2 x). A simple change of variable in the integrals leads to the equality 

4=Q| |/| | 2 , (B3) 

which proves the equivalence between the condition dBll and the wavelet admissibility condition l|6} on the plane. 

Second, let ty(ui) be a function in L 2 (S 2 ,dil). We show in a similar way that the wavelet admissibility condition (1131 on the 
sphere applied to is induced by the condition: 

0<4= / d X — / dSlo\(-*^ x .a\F)\ 2 (B4) 
Jo Jo a ~ Js 1 

for any F(oj) in L 2 (S 2 ,dn). Decomposing ^(uj) and F(uf) in spherical harmonics through JA4L we get the following relation: 
{^^xA F ) = E/ eN E|m|</ (*^„. x ,fl) /M /;m, with (* Wo , x . fl ) /m = E|„|</ D /™(^o, X)(^a)in (see relation ED). From the orthogonality 
relation \A5\ for the Wigner D-functions, the integral 1^ finally reads: 

ten \ m \<i 

This proves that the condition (IB4t implies the wavelet admissibility condition Jl 3i on the sphere. The condition iB4\ indeed re- 
quires that be strictly positive and bounded for all I by a strictly positive constant c, < < c, while the wavelet admissibility 
condition dl 31 reads < < oc. In the original group theoreti c approach, the condition dB4> defines the admissibility, through 
the square-integrability of the considered group representation ( Antoine & Vanderehevnst 1999). Our admissibility condition for 
the reconstruction of the signal is therefore slightly less restrictive. 

Establishment of the correspondence principle 

We can now prove the correspondence principle by showing that, if tp(x) in L 2 (M. 2 ,d 2 x) satisfies the reformulated wavelet 
adm issibility condition dBll . then the inverse stereographic projection ^(u) = [n _1 ^](w) in L 2 (S 2 ,dfl) satisfies the condition 
iB4\ (Boadanova et al. 2005), and therefore the admissibility condition Jl 3i . 

First, we establish the upper bound Inj-i^ < oo for any F(u>) ^0 in L 2 (S 2 ,dfl). The scalar product of functions in L 2 (S 2 ,d$V) 

is invariant under translation by u>o, rotation by x- Considering a unitary operator n between L 2 (S 2 ,dil) and L 2 (M. 2 ,d 2 x) implies 
by definition that the scalar product of functions is preserved by the projection. Moreover, for a unitary, radial, and conformal 
projection operator n, the conjugation relation J17l > holds between the dilation operator D(a) in L 2 (S 2 ,dfl) on the sphere and 
the dilation operator d(a) in L 2 (M. 2 ,d 2 x) on the plane. Indeed, the dilation operator d(a) of functions in L 2 (M. 2 ,d 2 x) on the plane 
defined by is also unitary, radial and conformal. The conjugate operator U~ l d(a)H on functions in L 2 (S 2 ,dil) on the sphere 
therefore also satisfies the same properties. Consequently, the uniqueness of the unitary, radial and conformal dilation operator 
in L 2 (S 2 ,d£Y) on the sphere (see appendix lAl implies that the operator n -1 ii(a)n is identified with the dilation D(a) defined in 
( I10i . We therefore get: 



/ ( n-'V)= / d P i at-^)(P) 

J 50(3) L 



(B6) 
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with 

f +oa da 

i(n-^)(P) = / — \{d(a)i/>\Wr l (p)F)\ 2 . (B7) 
Jo a 

Similarly, from the invariance of the scalar product of functi ons i n L 2 {R. 2 ,d 2 x) under translation by Xq and rotation by \ we may 
write the upper bound of the wavelet admissibility condition (IB 1 i on ip(x) on the plane as 

/■27T r+oo j p 

4=/ d W — / #%\{d(a)il>\r- 1 0dt- 1 (&)f)\ 2 <oo, (B8) 

JO JO fl JR 2 

for any f(x) ^ in L 2 (R 2 , d 2 x). By continuity of the integrand in the variables Xq and %, this finally implies that 

r+°° An 

\ —\(d(a)^\f)\ 2 <w, (B9) 
Jo fl 

for any f(x) ^0 in L 2 (RV 2 i). Consequently, for any F(uf) 4 in L 2 (S 2 ,dn), the function /(f) = [il/T^F ](i) is in L 2 (R 2 ,d 2 x) 
and differs from zero, and we readily obtain that i^j-i^ip) < oo, for any p e 5(9(3). The compactness of the group 5(9(3) finally 

ensures that Iqj-im < 00 f° r any f (w) 7^0 in L 2 (S 2 ,dfl). 

Notice that the choice of the measure of integration on scales da /a 3 is natural on the plane (see 0), while a priori arbitrary 
on the sphere (see J12t ~). It is however required to establish the correspondence principle between the wavelet formalisms on the 
plane and on the sphere, as it clearly appears from the former proof. 

Second, the lower bound < Iqx- 1 ^) ^ or an ^ F(uj) =^0 in L 2 (S 2 ,dtt) remains to be established. In that regard, we simply notice 
that the set of functions obtained by translation by u>q, rotation by x, and dilation by a of any non-identically null function is dense 
in L 2 (S 2 ,d£l). The lower bound of the wavelet admissibility condition on the plane iBl\ . < for any f(x) ^0 in L 2 (M. 2 ,d 2 x), 
ensures that a wavelet on the plane cannot be identically null. The unitarity of the stereographic projection therefore implies 
that the inverse stereographic projection ^(us) = [II" 1 ip](u>) is also different from zero in L 2 (S 2 ,dft). Consequently, the set of 
functions ^u) 0lXt a is dense in L 2 (S 2 ,dil), and considering any F(ui) ^0 in L 2 (S 2 ,dft), the scalar product u ,x,a\F) cannot be 
identically null in u>o, \, and a. Again, the continuity of this function in all its arguments ojq, x, and a ensures that it is non-zero 
on a set of non-zero measure in the corresponding Hilbert space. This strict positivity of the integrand in Iqj-i^ finally guarantees 

that < Iqj-im f° r any F(ui) ^ in L 2 (S 2 ,dST). This completes the proof of the correspondence principle. 

Uniqueness of the stereographic projection operator 

Let us prove that the stereographic projection defined in (1191 is the unique unitary, radial, and conformal diffeomorphism 
between the sphere and the plane. The projection operator II between functions G in L 2 (S 2 ,d£Y) on the sphere andg inL 2 (R},d 2 x) 
on the plane is generically expressed in terms the inverse of the corresponding projection operator ir between points on the sphere 
S 2 and on the plane M 2 , applied to the argument of the function considered. The projection tt is a radial diffeomorphism, that is a 
continuously differentiable bijection between S 2 and M 2 , which only relates the radial variables r on the plane and 9 on the sphere 
independently of tp, and which leaves ip invariant. In these terms, the operator II and its inverse II" 1 are given by the relations 
dl9t and i20\ . with tt{9, ip) = (r(8), ip), and its inverse ir~ l {r, p) = {0{r), <p). The functions r(9) and its inverse 9(r), together with 
the function fi(r), have to be fixed by the additional properties required. 

First, the projection tt is a mapping between the sphere and the plane, ir : u = (x l ,x 2 ) = (9, <p) £ S 2 — * x' = (x n ,x' 2 ) = (r, ip) £ M. 2 . 
It must be conformal in order to preserve the measure of angles and directions, which are defined locally in tangent plane at 
each point of S 2 and K 2 . This explicitly means that the metric g-,-(r, tp) = (dx* /dx")(dx l /dx'j)gki(9(r), ip) induced on the plane 

by the projection from the sphere is conformally equivalent to the Euclidean metric gij(r,ip) = diag(l,r 2 ) on the plane, with 
i,j,k,l £ {1,2}. The original metric on the sphere is the canonical metric gij(9,ip) = diag(l,sin 2 6>). The conformal equivalence 
reads by definition: g'ij(r,(p) = e^gij(r,ip), for some strictly positive conformal factor e^ r ', that is for a real function <f>(r). This 
condition requires that r{9) is a linear function of tan#/2. Let us recall that we consider the plane tangent to the sphere at the 
North pole (zo = 1). For the choice of coordinates (r, p) and (9, p) defined in §0] this gives uniquely: 

n 

r(0) = 2tan-. (B10) 

This radial diffeomorphism r(9) and its inverse 9(r) explicitly define the stereographic projection tt of points on the sphere S 2 
onto points on the plane R 2 , and its inverse. By stereographic projection, the sphere without its South pole is mapped onto the 
entire plane: r{9) : 9 £ [0, 7r[— > r £ [0, oof. A point u> = (9, p) on the unit sphere is projected onto a point x = (r, tp) on the tangent 
plane at the North pole, co-linear with ui and the South pole (see Fig. The conformal factor is also explicitly given by the 
condition of conformal mapping: e^ {r) / 2 = (1 + (r/2) 2 ) -1 , or equivalently e~^ {r{ - e ^l 2 = 1 +tan 2 (#/2). 

Second, the unitarity of the operator II on functions between L 2 (S 2 ,d£l) and L 2 (M. 2 ,d 2 x) identifies /i(r) with the conformal 
factor: fi(r) = e^ r \ This identity simply relies on the relation between the invariant Euclidean measure on the plane and the 
measure induced from the projection: r(9)dr(9)dp = e^MW sin9d9dip. 

In conclusion, the unique projection operator satisfying the properties required to ensure a correspondence principle between 
the formalisms of wavelets on the plane and on the sphere is therefore the stereographic projection operator dl9l l. 
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